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Abstract 

We perform a calculation on the polarization of J/ip production in deep 
inelastic scattering in the HERA energy range. For the inclusive production 
distributions, we find that the color-singlet contributions are consistent with 
the experimental data in the major region oi z {z > 0.4). Only in low z regions, 
there are some hints of the need of the color-octet contributions to describe 
the experimental data. For the polarization of J/ijj in DIS processes, we find 
the parameter a changes with Q^. Especially, at higher Q^, difference on a 
between color-singlet and color-octet contributions become more distinctive. 
In the two regions of lower and larger z, the polarization parameter a have 
different features. These properties can provide important information on the 
polarization mechanism for J production. 

PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 
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I. INTRODUCTION 



Studies of heavy quarkonium production in high energy colhsions provide important in- 
formation on both perturbative and nonperturbative QCD. In the conventional picture, the 
heavy quarkonium production is described in the color-singlet model (CSM) In this 
model, it is assumed that the heavy quark pair must be produced in a color-singlet state at 
short distance with the same angular-momentum quantum number as the charmonium which 
is eventually observed. However, with the recent Tevatron data on high px J / production, 
this color-singlet picture for heavy quarkonium production has become questionable. The 
observed cross section is larger than the theoretical prediction of the color-singlet model by 
a factor of about 30 ~ 50 0. This is called the J/ip {ip') surplus problem. On the theo- 
retical side, the naive color-singlet model has been supplanted by the nonrelativistic QCD 
(NRQCD) factorization formalism 0, which allows the infrared safe calculation of inclu- 
sive charmonium production and decay rates. In this approach, the production process is 
factorized into short and long distance parts, while the latter is associated with the nonper- 
turbative matrix elements of four-fermion operators. So, for heavy quarkonium production, 
the quark-antiquark pair does not need to be in the color-singlet state in the short distance 
production stage, which is at the scale of l/mg {mq is the heavy quark mass). At this stage, 
the color configuration other than the singlet is allowed for the heavy quark pair, such as 
color-octet. The later situation for heavy quarkonium production is called the color-octet 
mechanism. In this production mechanism, heavy quark-antiquark pair is produced at short 
distances in a color-octet state, and then hadronizes into a final state quarkonium (physical 
state) nonperturbatively. The color-octet terms in the gluon fragmentation to J/ilj{ip') have 
been considered to explain the J/iplip') surplus problems discovered by CDF . Tak- 
ing the nonperturbative {Oq^^{^Si)) and (Cg (^5'i)) as input parameters, the CDF surplus 
problem for J/ip and ip' can be explained as the contributions of color-octet terms due to 
gluon fragmentation. 

Apart from the NRQCD factorization approach (NRQCD FA) mentioned above, there are 
also other approaches for describing heavy quarkonium production in literature advocated 
in these years, such as the Color Evaporation Model 0, and the model of heavy quarkonium 
production with interactions with comoving fields [0]. 

Even though the color-octet mechanism has achieved some successes in describing the 
production and decay of heavy quarkonia, more tests of this mechanism are still needed. 
Recently, the photoproduction data from HERA put a question about the color-octet 
predictions for the inelastic photoproduction of J/ ip [0)0] (possible solutions for this prob- 
lem have been suggested in |T2|-p!5|). Most recently, the CDF collaboration have reported 
their preliminary measurements on the polarizations of the promptly produced charmonium 
and bottomonium states |]TB|, which appear not to support the color-octet predictions that 



the directly produced S'-wave quarkonia have transverse polarization at large pr jr^Jlg] . This 
discrepancy between the experimental results and the theoretical predictions shows that the 
predictions of NRQCD color-octet mechanism on the polarizations of heavy quarkonium 
productions may be questionable. However, the final conclusion about this problem can 
be achieved only if the polarizations of charmonium production in other processes are also 
measured and compared with theoretical predictions. In this paper, we will study the po- 
larization of J/ip production in deep inelastic scattering (DIS) processes at HERA collider. 
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The relevant photoproduction processes have been studied in |T^, and the leading order 
color-octet J/ip production and polarization in DIS have been studied in |20| (from 2 — >■ 1 
virtual photon subprocesses). These processes contribute to J/ip production in the forward 
direction. In this paper, we will complete these studies by calculating the J /ifj production in 
DIS from the NLO color-octet processes, i.e., from the virtual photon 2^2 subprocesses. 
With these calculations, we can study the z distributions of J / ip production and polarization 
in DIS processes at HERA, which will compensate the previous studies in the photoproduc- 
tion process. Moreover, since the photon virtuality can be large, electroproduction is a 
better process from which to test the color-octet mechanism and to extract the NRQCD long 
distance matrix elements than photoproduction. The latter process lacks any large energy 
scale other than the charm quark mass, and consequently, higher order perturbative correc- 
tions to leading order calculations are expected to be large. In addition, nonperturbative 
effects, such as higher twist corrections to the parton model, are less effectively suppressed 
in photoproduction than in electroproduction at large Q^. 

The rest of the paper is organized as follows. In Sec. II, we will give the polarized cross 
section formula for the inelastic J/tp production at the electron-proton collider. Here, we 
adopt the Weizsacker- Williams approximation to calculate the electroproduction cross sec- 
tion with the photoproduction cross section. The numerical results are given in Sec. HI. We 
will display the polarized cross sections for both the photon-proton collisions and electron- 
proton collisions. In Sec. IV, we give the conclusions. 



II. POLARIZED CROSS SECTION FORMULAS 

The polar angular distribution in J/ip decay is given by 

dV 



dcosO 



oc 1 + a cos 6, (1) 



where 6 is the angle between the lepton three- momentum in the J/ip rest frame and the 
polarization axis, a is the polar angle asymmetry parameter, 

where 

^ dajep ^ eJ/jjjX = 0)X) 
Exd(r{ep^eJ/i:{X)X)- 

Here A is the helicity of the produced J/ip. A = means J/ip is longitudinally polarized, 
and A = ±1 transversely polarized. 

In the Born approximation, the electroproduction cross section cr(ej9 eJ/ipX) is related 
to the 7*p cross section by 

Here Tj^l are the flux factors of the transversely and longitudinally polarized virtual photons 
respectively. 
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_ a(l + {l-yf) _ ail-y) , 

where y = k^* ■ P/kg ■ P is the fraction of the lepton's energy lost in the proton rest frame. 
In the typical kinematic region of HERA experiments, the difference between Ft and Fj^ is 
negligible. So, in practice, we can simplify the cross section formula of Eq.(Q) as 

dajep eJ/ipX) 

^y^Q2 = TTatot{l P ^ J/VX), and, atot = (Tt + (^l- (6) 

(Tt and are the cross sections for the transversely and longitudinally polarized virtual 
photons respectively. 

There are two types of contributions to the 7*p cross section: the direct photoproduction 
and the resolved photon production. The latter contribution is through the partonic content 
of 7* in the reactions. The cross section formulas (polarized and unpolarized) for the relevant 
partonic processes of these resolved photon production processes can be found in [|l9l . In the 
following we will calculate the polarized cross sections for the direct virtual photon {Q"^ > 0) 
processes, which include the following partonic channels, 

^^ + g^ic-c)['s['\'s['\'s!,'\'PP]+g; (7) 
7* + q/q - icc)['S['\'Si'\ ^'Pf ] + q/q. (8) 

In this paper, we calculate the above 2 —>■ 2 subprocess contributions to production of J/ip, 
and present the z distribution of the production rate. For this purpose, in the following we 
will not consider the 2 — > 1 subprocess contributions, because these contributions only take 
place in the forward region. 

To calculate these virtual photon subprocesses, we employ the helicity amplitude method. 



Following we choose the polarization vectors for the incident and the outgoing gluons 

as 

±) 



2 - Are[^2 ^3 #(1T75)+ #^3 ^2(1±75)], (9) 
= iV,[^3 # ^2(1 T 75)+ ^2 ^ ^3(1 ± 75)]. (10) 

Where q = Pi + 2^-j^P2 , Pi = —Q^ and pi is the momentum for the incident photon. p2, Ps 
and 62, 63 are the momenta and the polarization vectors for the incident gluon and outgoing 
gluon respectively. The normalization factor Ng is 

= , ^ (11) 

y'2(Q2M2 + su)t 

where M is the mass of J/ip. The Mandelstam invariants s, t, u are defined as 

s = (Pl+P2)^ i={p2-p3y, u={pi-p3y, (12) 
and they satisfy the relation 

For the transversely polarized photons, we have the following projection operator. 
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= <(T)e.(T) = -g,^ + P^'^^^ + P^^'^^ , (13) 

P2 ■ Q 



and for the longitudinally polarized photons, we have 

Pl ■P2 Pi- P2 



P;!; = e;{L)e,{L) = — (pi, + ^P2,){piu + z^P2u). (14) 



Providing with the above defined polarized vectors and polarized projection tensors, we 
can calculate the amplitude squared for the partonic processes with the definite helicities 
of the incident and the outgoing partons. After summing up the helicities and colors, we 
express the amplitude squared as the following form 

n 

where the short-distance coefficients F can be written as 

+Bi,[n][e*{X)-pie{X)-pi] 
+Q,[n][e*{X)-p2e{X)-p2] 

+D,,[nWiX)■p^MX).p,]. (16) 

Here n denote the intermediated states, which include ^S[^\ ^S[^\ ^Sq^\ ^Pf^- The 
A, B, C, D functions for different partonic processes are listed in the Appendix. 

For convenience, we consider the J/ip polarization in the target iraiae. In this frame, the 
covariant expression for the polarization vector of J/ ip with helicity A = reads [O] 

^r-"'M-}iv..-^p.. (IT) 



With this equation and Eqs.([T^|T6|), we can calculate the polarized cross sections, and obtain 
the polar angle asymmetry parameter a in Eq.(|^). 

We have checked that our cross section formulas can reproduce the photoproduction cross 
section formulas [jl9| at = 0, and that the unpolarized cross section for the color-singlet 
process is consistent with that of 



III. NUMERICAL RESULTS 

The production rate of J / ip and its polarization parameter a depend on the sizes of the 
NRQCD long distance matrix elements, including the color-singlet matrix elements and the 
color-octet matrix elements. We choose these matrix elements as 



{Oti^Si)) = l.WGeV^ (18) 

(Oti^Si)) = 1.06 X 10-^GeV\ (19) 

{OtCSo)) = 3.0 X lO-^GeV^^^ (20) 

{OtCPo))/ml = 1.0 X lO'^GeV^ (21) 

{OtePj)) = i2J + l){OtePo)). (22) 
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The color-singlet matrix element {Ofi^Si)) can be related to the cc wave function at the 
origin, and can be taken from the leptonic decay width of J/ip. The last equation comes from 
the approximation of heavy quark spin symmetry of NRQCD. The value of the color-octet 
matrix element {Og{^Si)) is taken from a fit to the large pt J/ip production at the Tevatron 



25| . This matrix element is not important to J/ ip photoproduction both for the production 
rate and the polarization parameter a. On the other hand, the other two color-octet matrix 
elements, {Ot^So)) and (Oti^Po)), are known to be very important in the inelastic J/ijj 
photoproduction ||l0| , p!9| . However, their values are not well determined from the present 
experimental data on J/ip productions. Here, we just follow Ref. []TU[ and take their values 



tentatively as listed above (which are also consistent with the naive NRQCD velocity scaling 
rules) to see what are their contributions to the cross section and the polarization a of J/ijj 
production in DIS at HERA. For the numerical evaluation, we choose rric = 1.5 GeV, and 
set the renormalization scale and the factorization scale both equal to /x^ = (2mc)^ -|- Q^. 
For the parton distribution functions of the proton, we use the Gliick-Reya-Vogt (GRV) LO 
parameterization . 

We first display the z distribution of the inelastic J/ip production rate in DIS region 
at HERA, comparing the theoretical predictions with the ZEUS data Q. The kinematic 
region is 2GeV^ < < SOGeV^ and AOGeV < W^,p < ISOGeV. In Fig. 1, the dotted 
line is for the color-singlet contribution, the dotted-dashed line for the direct virtual photon 
contributions from the color-octet processes, and the dashed line for the resolved virtual 
photon contributions from the color-octet processes. For the above two components of 
the color-octet contributions, we choose the color-octet matrix elements of {Of{^So)) and 
(C^(3Po)) as {Oti^So)) = {Oti^Po))/ml = O.OOSGeV^ And the sohd lines correspond to 
the NRQCD FA predictions (including the color-singlet contributions and the color-octet 
contributions) for the two choices of the color-octet matrix elements: (I) for the lower solid 
line, 

{OtCSo)) = 3.0 X 10-^GeV\ 
{Ot{'Po))/ml = 0; (23) 

(II) for the upper solid line, 

{Oti'S,)) = 0, 

{OtCPo))/ml = 1.0 X IQ-^GeV^. (24) 

^From this figure, we can also see that the rapid increase of the color-octet predictions at 
large z is not supported by the experimental data. This is similar to the case of the J/i/j 
photoproduction {Q^ = 0), and may indicate that at large z the calculations of the color-octet 
processes are also unreliable due to higher order v"^ contributions |T^. On the other hand. 



we find that the color-singlet contributions are consistent with the experimental data in the 



^These data have also been compared with the theoretical predictions in |26|. However, this study 
is not relevant to the electroproduction process because the author made an approximation for the 
cross sections of virtual photon processes by using those of real photon processes. 



6 



major region of z. However, as indicated by the color-singlet photoproduction processes 



calculations |TT|,[19[, the color-singlet cross sections have large theoretical uncertainties due 
to input parameters such as the normalization and factorization scales and the charm quark 
mass. In this context, it is difficult to immediately conclude that the experimental data 
are just saturated by the color-singlet contributions and there is no more space for the 
color-octet contributions. Especially, in the lower z region the experimental data lie above 
the color-singlet contributions, where the dominant contributions come from the color-octet 
channels in the resolved photon processes. (The color-singlet contributions in the resolved 
photon processes is much smaller than the color-octet contributions, and is not presented 
in Fig. 1). This may be viewed as a hint of the need of the color-octet contributions in 
describing the experimental data on J/?/; photoproduction in the lower z region. 

We now turn to study the angular distributions of J/?/' productions. Since the decay 
angular distribution parameter a is normalized, its dependence on parameters that affect 
the absolute normalization of the production cross sections, such as the charm quark mass, 
the strong coupling constant, the renormalization and factorization scales and the parton 
distribution functions, cancels to a large extent and does not constitute a significant un- 
certainty. That is to say, in some sense, the polarization parameter as a tool for testing 
production mechanism is more efficient than the absolute production rate. In Fig. 2, we 
first display the polarization parameter a as a function of z for the 7*p processes, where the 
cm. energy of 7*p system is set to be W^*^ = lOOGeV. Here in this figure we do not include 
the resolved photon processes contributions. The solid lines are for the CSM predictions, 
and the other two lines are for the NRQCD FA (including both the color-singlet and color- 
octet contributions) corresponding to the two choices of the color-octet matrix elements: the 
dashed lines for the choice of Eq. (p3D and the dotted-dashed lines for the choice of Eq. (p4|). 
In order to see the dependence of the polarization, we choose four typical values for Q^: 
0, 4, 10, AOGeV^. From these plots, we can see that the polarization parameter a in large 
z region change from positive values to negative values as increases. This means that 
in large z region J/ip production will be dominantly longitudinally polarized at high Q^. 
Especially, the polarization parameter a for CSM predictions are more sensitive to than 
those for the NRQCD FA predictions. At high Q^, e.g., = AOGeV'^, the difference on a 
between CSM and NRQCD FA can be distinguished for large z J/ip production (see Fig. 2). 
In this case, the CSM prediction approaches to a = —0.7 while the NRQCD FA predictions 
are both above —0.4 for the two choices of the color-octet matrix elements (Eqs. ( p3D and 
(1^)). In addition, we note that if ^Sq^^ channel dominates J/ip production the polarization 
parameter a will be close to zero because in this case the produced J/ip are unpolarized, 
which can also be seen from Fig. 2. 

In Fig. 3, we show the parameter a as a function of z for ep collisions at HERA (including 
the resolved photon contributions). As in Fig. 2, we choose two typical regions for Q^: (a) 
2Ge1/2 < g2 < SOGeV^, and (b) lOGeV^ < < 80GeV\ From this figure, we also 
find that the J/ip polarization changes with Q^. Especially, at higher Q^, the difference 
on a between CSM predictions and NRQCD FA predictions become more distinctive. In 
the two regions of lower and larger z, the polarization parameter a have different features. 
In the lower z region, the NRQCD FA predicts J/ip being transversely polarized for the 
octet matrix elements choice of Eq. (p4|), while the CSM predicts J/ip to be only slightly 



polarized (almost compatible with unpolarized). We note that in this region, the dominant 
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contributions in the NRQCD FA come from the color-octet resolved photon processes (see 
Fig. 1), which have the similar properties as they have for J /ip hadroproduction, contributing 
to J /ip transverse polarization In Fig. 3, for larger z, on the other hand we find that the 



CSM predicts J/i}) being more polarized, which is however longitudinal. The NRQCD FA 
predicts J/^ less polarized (also longitudinally if ^Pf^ dominates). Again, we note if J /ifj 
production is dominated by ^Sq''' channel, it will be unpolarized as that for YP processes. 

In the above analyzes, we present the NRQCD FA predictions of the polarization param- 
eters a by using the color-octet matrix elements as in Eqs. (^) and (p^ . However, there 
are other sets of parameterizations for the two color-octet matrix elements {0^{}Sq)) and 
(Og (^Pq)) in the literature |1^,|13|. With these rather small values for the two color-octet 
matrix elements the above predictions of a in the NRQCD FA will be changed. For 

example, if we follow the values obtained in [1^, the two equations of ( ^3]) and (^4]) will be 
changed to 



for case (I), and 



{Oti^S^)) = 5.72 X IQ-'^GeV^, 
{OtePo))/ml = 0; (25) 



{OtCSo)) = 0, 

{OtCPo))/ml = 1.62 X 10-^GeV\ (26) 



for case (II). If we use the above values for the matrix elements of (OgC^So)) and {Og{^Po)), 
the difference of the polarization parameter a between the CSM and the NRQCD FA will 
be reduced, because with these rather small values for the color-octet matrix elements the 
color-octet contributions to J ftp production are much less important than those with pa- 
rameterizations of Eqs. (^) and (p^). This influence is presented in Fig. 4, where we show 
the same plot as in Fig. 3 but with the new parameterizations of the color-octet matrix 
elements as in Eqs. (^) and (PB|). From this figure, we can see that the difference between 
these two approaches is reduced, though in some regions there are still some differences. 



IV. CONCLUSIONS 

In this paper, we have calculated J/ip production and polarization in DIS at ep colliders 
in the energy region relevant to HERA. Compensating for the previous color-octet leading 
order calculations, we have calculated the cross sections for the 7*p 2 — >• 2 subprocesses 
which are needed for theoretical evaluations of the z distributions of the production rate 
and the polarization parameter a. For the inclusive production distributions, we find that 
the color-singlet contributions are consistent with the experimental data in the major region 
oi z {z > 0.4). Only in the low z regions, there are some hints of the need of the color-octet 
contributions to describe the experimental data. 

For the polarization of J/ip in DIS processes, we find the parameter a changes with Q^. 
Especially, at higher Q^, difference on a between the CSM and the NRQCD FA become 
more distinctive. In the two regions of lower and larger z, the polarization parameter a have 
different features. In the lower z region, the NRQCD FA predicts J/ip being transversely 
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polarized if the ^Pj channel dominates, while the CSM J/ ip being almost unpolarized. On 
the other hand, in the larger z region the CSM predicts J/ ip being longitudinally polarized, 
while the NRQCD FA predicts J/ip being unpolarized if the -^Sq^^ channel dominates. In 
conclusion, the polarization of J/ ip in DIS processes at ep colliders will give another indepen- 
dent test for the production mechanism and help to clarify the present problem concerning 
the disagreement of the theoretical predictions with the experimental data on polarization 
in hadropro duct ion processes at the Tevatron. 
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APPENDIX: 

In this appendix, we list the A, B, C, D functions of Eq.([T6|) in the short distance 
coefficients of the amplitude squared for every partonic processes^. For convenience, we 
define the following variables relevant to the Mandelstam invariants (s, t, and u): s = 
2pi-p2 = s + Q^, u = -2pi ■p3 = u + Q'^, t = 2p2-p3 = i 
^* + ^_cc(35i,l) + (7: 

16M(47r)3aa^eg(0f(3g,)) 

a = s'^[stu{s + t + u) - {st + tu + suf] - 2Q^{Q'^t'^ + stu){s'^ + t^), (A2) 

b = 2s'^[2QH^ - (s^ + t^){s + t + u)], (A3) 

c = 2[4QH^ + 4QH{su -t^ - tu) + 2Q'^s{st^ + sv^ - sH - t\ - tv^) 

-s\s + t + u){s^ + u% (A4) 

d = As[QH{2QH -s^ -t^ - tu) - s^{s + t + u)\. (A5) 



^It is easy to check that our A, B, C, D functions reproduce the results of Ref. |jT^ in the 
photoproduction limit. Using these expressions, one can also get the unpolarized cross sections for 
different subprocesses, which were reproduced by the authors of for the 7^ + g{q) processes. 
However, there is a factor of 2 difference from theirs for the "ii+giq) processes except the 72 + 5 ^ 
cc(^S'o,8) + g process. We note, however, that our unpolarized cross sections for the color-singlet 



processes 7j'/^ + 5 — ^ cc( 5*1, \) + g are consistent with those of ^2|. The computer code for these 
expressions can be obtained by requiring to ktchao@pku.edu.cn. 
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72 + (? ^ cc{^Si, 1) + g: F is the same as Eg. {\Kl\) and 6 = 0, 

a = -2QH%QH + su), (A6) 
c = 4Q2[2Q2(qV ^ _ ^3 _ ^2^^) ^ ^^(^^ - t)(t + m)], (A7) 

d = 4Q^st[su + t{Q^ - M% (A8) 

7tl + cc(^S'i, 8) + g: the functions are the same as the above 'Jj- l~'<~ 9 ^ cc(^Si, I) + g 
processes but muhiphed by the factor 

8 (Ofi^Si))' ^ ' 

1t + 9 ^ cc(^S'o,8) +g: 

MsH[s + tY[s + 'u)^(t + -u)^ 

+s\s'' + t^ + u^ + {s + t + u)% (AlO) 
72 + ^ cc(i5o, 8) + 

+s2((t^ + M^)(t2 + 2tM + 2u^) + 2s(t + u){t^ + tu + u^) + s2(t + uf)]. (All) 
7^ + cc(3Pj, 8) + g: 



MHH\s + tf{s + uY{t + uf' 
a = -t[-8Qh^{s^ + t^){s^ - tu){st - u^) + 2QH{s^{hr' + 3t^M - 8tV + IQt^u^ + 2m^ 



(A12) 



-sH{f + 9t^M + 4t2M2 - lOtM^ - 2u^) + /(7t^ - 7t^M + 4t2M2 - 2tu^ + 4m^) - st^u{<ot^ 
+ llt^u - fu^ + htv!^ - u^) + s5(7t=^ + - 2tM2 + 8^3) + 45^(2^^ + ^x^) + 2(t + m)s^ 
-t\^ + 3t^u^ + 5t^M^ + t^v!') - 2Q\s\ht^ + It'u - 5t\^ + 6t\^ - 8t\^ + htu^ - 2u^) 
-sh'^if + Sfu + m'^u^ + 5t^M^ + lOt^M^ - 5tM^ - M*') + s\9r> + 2t^M + 4^3^^ - St^M^ 
+8tM^ - 4m^) + sHu{t^ + 7t^M + t^u'^ + 7t\3 - tM^ + u^) + s^(9t^ + 7t\ - 2t^u'^ + 8tti^ 
-6n^) - (t + u)st^M(2t^ + 8t^M + ht^u^ + 9tu^ + 2u^) + s^(9t^ + 4^^^ + htu^ - 4^^) 
+s®(3t - 2m) (t + m) - tV - 2tV - tV) + g2s(s^(3t^ + 29t^M + 56tV + 34t^M^ 
+26^3^^ + Qt^u^ - 6tu^ - 6m^) + s^(llt^ + 30t^M + Gt^u'^ - 8t\^ - 22fu'^ - 13tM^ - 18m^) 
+sHu{St^ + 45t^u + 76t\^ + 60t^u^ + AQt^u* + lltu^ - u^) + 2s^(7t^ + U\ - ?,t^u^ 
-Ut^u^ - lltu^ - Uu^) + 5*^(14^^ + 13t^M - 8t^M^ - 19tM^ - 28m^) + st^u^it + M)(7t^ 
+2At^u + ISt^u^ + 18tM3 + 5m^) + s^(llt^ + t^u - I2tu^ - ISv!^) + 2{t + m)^^^^^^ 
+3(t + u){t- 2u)s^) + /^^(lOt^ + ASr'u + 80t^u^ + 67t^u^ + Aht^u^ + 19tu^ + 3u^) 
+s^u(6t^ + 48t^u + im%^ + 114t3M3 + 85^2^^ + Ahtu^ + 12^^) + s^{t^ + 22t^u + 80t^M2 
+ 114^3^3 + i02t\4 + 61tM^ + 22u^) + s^(2t^ + 27t^M + Qlt^u^ + 85^2^3 + 61tM^ + 26m^) 
+s^tM^(t + u){<6t^ + 16t^M + llt^M^ + 8tu^ + 3u^) + s^{t^ + 19t^u + 45t\^ + 45tu^ 
+22u^) + s\(llt2 + 19tu + 12^2) + s2tV(t + m)2 + s^°u3(t + m)] (A13) 
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b = 2[4QV(st - u'^){2sH - s\ - 2shu - s^u^ - AsH\ - 2sHu^ - st\ - 2st^u^ - t^u'^) 
-2Q^st{s^{2t^ - 2t% - ISt^u^ - let'^u^ - 5tu^ - 2u^) + sH{f - t^u - ISt^u'^ - I2fu^ 
-lOtM^ - 2u^) + s^(8t^ + 9t^M - St^u^ - 3tu^ - Au^) + s^(llt^ + Gt'^u + tu^ - Au^) 
+st\{t + u){t^ - At^u - tv^ - 2u^) + s^{Qt^ + 5tM - 2u'^) - 2t^u^{t^ + tu + v^) + 2sH) 
-sH{f + IQt'u + 25t^u^ + 23t^u^ + 25t^u^ + 16tu^ + 4n^) - sh{f + lOt^u + 2lt^u^ 
+m'^u^ + 17tM^ + I2u^) + s^(4t^ - t^M + 2fu^ - 3tu^ + 2u^) - sH\{t + u){2t^ + 9t^u 
+3fu'^ + 8tu^ + 4m^) + sH{2t'^ - ht^u - l^t^u^ - 3tu^ - IOm^) 

+s^{3t^ + t^u - 2tu^ + 2u^) - sh^u^if + 2tu - u^) + s^t{t - u)]M'^ (A14) 

c = -2[-16Q*^t^(s^ - tn)(st - u^) - AQH'^{s^{t^ + 3t^u + Atv!^ - 4m^) + su{2t - u)(3t^ 
+5t'^u + 2tu^ + 2u^) - s^t + 2m) (St^ - + 4^2) - As^{t^ + u^) - 2{t + + tV 
-5t\^ - 6t\^ - 2tu^) + AQH{s'^{t^ + 8t^M + 9t\^ + 8t^M^ + 10t\^ + tu^ + 2m*^) 
+stM(2t^ + m% + llt^M^ + I3t^u^ + 7tM^ + u^) + s^uit'^ + t^u + Qt^u^ - 3tM^ + 4u^) 
- t^u - St^u^ + 2tu^ - 6m^) - s^{2t^ - 3t^u + 2tu2 - Au^) + s^(t^ + + 2u^) 

+t^U^ + 2tV - t^U^ - t^U^) + 2Q2st(s^^3(^^^^4 ^ ^Q^3^ ^ g^2^2 _ ^^3 _ ^u^) + sV(llt^ 
+9t^M + t^M^ + StM^ - hu^) + S^(5t^ + ht^U - 12t\2 _ g^^S _ ^§^4^, ^ ^3^(5^4 _ ^g^2^2 

-5tM^ - 12u^) + 2s^(2t^ - t^u - tu^ - 7u^) + ^'^(St^ + 3tu - 5u^) - 2tu^{2t^ + 2tu + u^) 
+2s^t) + s^u{2t' + 18t^M + 6t\2 - ISt^M^ - tM^ - 2u^) + s*^M(13t^ + I2t^u - Iht^u^ 
-6tu^ - 6u^) + sHu^iAf^ + ISf'u + 12t\2 - tu^ + 3u^) + sHu^{2t^ + 13t\ + 7t\^ 
+tu^ + u^) + s\t - u){t^ + 8t\ + 7tu^ + 6u^) - s^{2t'^ + tu- 2u^){t - u) 
+sh'^u\t^ - 2tu - u^) - (t - u)sh]M^ (A15) 

d = -2[-8Qh'\st - M^)(2s^ - s\ - 3stu - su^ - tu^) + AQH'^ {stu^ {t^ + 2t^u + Qtu^ + u^) 
+sH{bt^ + 7tu + 4^2) - s^uit + M)(2t2 + 3tu + Sm^) - s2M(2t2 + u^){f + 4tM + u^) 
+s^(5t + 2u)(t + n) + 2s'^(2t + u) + tV + t^u^ + 2tV) + 2Q2st(s2(2t^ + lAt^u + STt^w^ 
+42t\=^ + 37t\^ + IStM^ + 2m^) + stu{At' + ISt^w + 3lt^u^ + 33t\^ + 25tM^ + 6m^) 
+s3(3t^ + 21t^M + 3Qt^u^ + 22t\3 + mu^ + 3m^) + s^(7t^ + m^u + St^^^ + Utu^ - u^) 
+s^{7t^ + m'^u + lltn^ - 3n^) + s^{7t'^ + I2tu - u^) + 2t^u^{t^ + tu + 2u^){t^ + tu + u^) 
+AsH) - s^if - 7t% + 3t^u^ + 29t\^ + 6tu^ + Au^) - s^(3t^ + At^u + 19t\^ + lOtu^ 
+8u^) + sH'^u^{t^ + At^u + l^t^u^ + mu^ + 6m^) - sH{t + u){t^ - t^u - Iht^u^ + Otw^ 
-8m^) -{t + u)sHu{t'^ - 5t^u - 15t\^ - 12tu^ - 5u^) - s\5t^ + At^u - tu^ + Au^) 
+sH^u^{f + t^u + Stu^ + M^) - 2{t - u)s^t]M^ (A16) 

+ cci^Pj, 8)+g: F is the same as Eq.( |A12| ), 

a = -QH^[-8QH^{s^ - tu){st - u^) - 2QH'^{su{6t'^ + llt\ - t^u^ + htu^ - u^) + ^^(t^ + 9t^u 

+t\^ - 3t\^ - 5t\'^ - tu^) + 2Qh{s^{t^ + U^u + llf^u^ + 6t^u^ + 16t\^ + 4tM^ + 2u^) 
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+8t\ + 5t\,^ + 9t?/ + 2u'^)stu(t + u)- s\8t^ + 7t\ + 6tu^ - 5u^) - {3t - u){t + u)s^ 
+fu^ + 2t^u^ + iV) + s{s + u){s^{t + u){3t'^ + u^) + s^{9t^ + Ut^u + m'^u'^ + 2tu^ + 2,u^) 
+s^(9i^ + 2At% + 2>2t\'^ + 20^^1*2 + 7tu^ + Au^) + s'^{2>t^ + m^u + AAt'^u^ + 42t^ii^ 
+27t2M4 + 8tM^ + 2^6) + Ssi^-u + 2^st^u^ + 35siV + 25siV + 125^2-1^5 ^ 25^1^^ 
+2iV + 4iV + 2iV)], (A17) 

6 = -4(s + tf{s + u)t'^M^Q^[2QHu{st - u^) + s{sH'^ - s^u^ + st^ - stu^ - 2su^ - 2t^u^ 
-2tu^ - 2u% (A18) 

c = -4Q2tM2[-8Q^t=^(s2 - tu){st - u^) - 2QH{s^{t'^ + 3t^u + 4t?/ - Au^) + (2t - u)su{'it^ 
+bt^u + 2^^^ + 2u^) -{t + 2u)s^{5f -tu + Au"^) - As'^it^ + u^) - 2{t + u)s^ + t%'^ 
-5fu^ - Qt^u^ - 2tu^) + 2Q'^{s^{f + St'u + Qt^w^ + m^u^ + 9t\^ + tu^ + 2m^) 
+stM(2t^ + m\ + llt^u^ + Ut^u^ + 7tM^ + -u^) - s'^if - t^u - At^u'^ - tu^ - Qu*) 
-s^u{t + 2m) (t^ - St^M + 21%^ - 2^3) + s'>u{t^ + to, + 4?7,2) + 2s^u{t + u) + t^i? + 2t^u^ 
-t^u^ - t\^) - s{s + u)(3sh'^ + 6s^t?i + 3s^u^ + AsH^ + 9sH\ + lOs^to^ + 5sSr^ 
+6sH^u + ISs^t^M^ + 14s-'^to''^ + 7s^u^ - s^t^ + 2sH^u + 12sH\'^ + ISs^t^?/ + 5s^to^ 
+352^^ - st^-u - Gst'^u^ - 6st^u^ - Ast'^u'^ + stu^ + 2su^ + At^u^ + Afu^ + 2tu^)\ 

d = 4Q2t[4Q4i2(^st - u^){2s^ - s^u - 3sto - sm^ _ t^^) + 2QH{{2t* + 8t^M + t^w^ ^ 4^^3 

+M^)s^M - s^(5f'^ + 7t^u + 4m^) + s^u{Af^ + 5t^?i + 6tM^ + u^) - sto^(t3 + 2t^u + 6tM^ + -u^) 
-s\7t^ + 3t?i + 4n2) - 2s\t + u) - - t^u^ - 2t\^) + s{s + u){s\3t'^ + 2tu - u^) 
+s^{6t^ + At^u + Atu^ - 2u^) + s3(t^ + 2t^u^ - 3m^) - s'^{At^ + I2t% + 12t3M2 ^ 14^2^3 
+8tM^ + 2m^) - 2st^ - Ust^u - 21stV - 24st^'u^ - 19st\^ - 6stu^ - 2fu - At^u^ 
-Sfu'' - 6t^u* - At^')] (A19) 

j^ + q^cc{^Si,8) + q: 

{Anfaay^iOteS,)) 

a = {Qh + su){2Q^{s + tf - 2Q^s[{2t + s){s + t + u) - tu] + s'^[{s + tf + {t + uf]}, (A21) 

b^As\Q^-u)[Q\s + t)-sM^], (A22) 

c = 4[-2g^t + 2Q^{s^ + 3st + f + tu) - AQ^st{t + u + s)+ Q'^s^{2f - - - Asu) 
+2{t + u)s^u + 2/^] , (A23) 

d^8s''{Q'' -u){Q^ + QH-sM^), ( A24) 
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where Cq is the electric charge of the hght quark q. 

ll + Q^ cc{^Si, 8) + q: F is the same as Eg . (|A20|) and b = 0, d = 0, 

a = 2QH{s + tf{Q'^ - sf, (A25) 

c = SQhM^{Q'^ ~ sf. (A26) 
7^ + g ^ cci^So, 8) + g: b = c = d = Q 

72 + g cc{^So, 8) +q: b = c = d = 

7^ + g ^ cc(3Pj, 8) + g: 

16(4vr)3aageg(0s^(3Po)) 

a = t[-8(s2 + t2)Q6^ + 2Q^(s2(5t2 + Atu + 2u^) + 2st^{2t - u) + 2s^ + 4s^t + + 4t='u 
+tV) + 2Q'^s{su{t - u){M + 3m) - s^(3t + 5m) + tu{2t + m)^ - s^u{t + 5m) - 3s^) 
+s2(M^(4t2 + 6tM + 3m2) + 2sM(4t + 5m) (t + m) + 2s^(3t + 5m) 

+2s^(t + m) (2t + 7m) + 3s^)] (A30) 
6 = -8M^[Q2t(s2(3t + m) + stit + m) + + t^n) - s^(s + m)((s + t)^ + tM)] (A31) 

c = %M\-^Qh^ + 2g^t((s - t)2 + (t + m)2) + g2t(4t5^ _ (s + _ ^)2) 

+s2(M + s + 2t)(M + s)^] (A32) 

(i = %M\-2QH^{?,s + m) + g2t(sy(4t + m) - 2(t - u)s^ + + 2tM2) 

+s2(s + m)(2s(m + s) +t(M + 3s))] (A33) 

7l + 9 ^ cc(3Pj, 8) + g: is the same as Eq. ([A29| ) 

a = 2Q¥[-4QV + Qh{2s{2t - u) - 3s^ + {2t + uf) 

+s{s + m) (4st + SU + (2t + m)^)] (A34) 

6 = -8QVM2(s + M)(s + t) (A35) 

c = 8Q2tM2[-4Q^t + 2Q2((s - t)^ + (t + m)^) - (s + m)((s - m)^ - Ast)] (A36) 

d = -8QHM^ [2QH{3s + m) - (s + 2t) (s + m)^] (A37) 
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Figure Captions 



FIG. 1. The differential cross section da/dz as a function of z for J/ip production in DIS at 
HERA: 2GeV'^ < < SOGeV^, AOGeV < W^,p < ISOGe^. Tlie dotted line is for the color- 
singlet contributions, the dotted-dashed line for the direct virtual photon contributions from 
the color-octet processes, and the dashed line for the resolved virtual photon contributions 
from the color-octet processes, where octet matrix elements take values as (Og (^Sq)) = 
{Ot{^Po))/ml = O.OOSGeV^. The sohd lines correspond to the total cross sections for the 
two choices of the color-octet matrix elements: (I) the lower solid line for choice of Eq. (^3|), 
and (II) the upper solid line for Eq. (|2^). 

FIG. 2. The polarization parameter a as a function of z in 7*p processes, where W^p = 
lOOGeV^ with four typical values for Q^. The solid lines are for the GSM predictions, and 
the other two lines are for the NRQGD FA predictions (including both the color-singlet and 
color-octet contributions): the dashed lines for the choice of Eq. ( PB| ) and the dotted-dashed 
lines for the choice of Eq. (^) . 

FIG. 3. a as a function of z for J/ip production in DIS at HERA. The definitions of the 
curves are the same as those of Fig. 2. 

FIG. 4. The same as the plot of Fig. 3, but with different parameterizations for the color- 
octet matrix elements (here Eqs. (pSj) and ([26| ) are used). The solid lines is the GSM 
prediction, and the other two lines are for the NRQGD FA predictions: the dashed line 
for the parameterization of Eq. (|25| ) and the dotted-dashed line for the parameterization of 
Eq. (ID. 
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